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Abstract. This article concerns the stabihty of a model for mass-spring sys- 
tems with positive damping and negative stiffness. It is well known that when 
the coefficients are frozen in time the system is unstable. Here we find condi- 
tions on the variable coefficients to prove stability. In particular, we disprove 
the believe that if the eigenvalues of the system change slowly in time the sys- 
tem remains unstable. We extend some of our results for nonlinear systems. 



1. Introduction 

In this article, we present conditions for the stabihty of solutions to the differ- 
ential equation 

u"{t) + b{t)u'{t) + k{t)u{t) ^ 0, (1.1) 
■w(io)=uo, u'(to)=wi, (1-2) 

where the coefficient k may have negative values. This equation has been used 
for modelling mass-spring systems, where the mass is one unit, the coefhcient b 
produces a damping effect proportional to velocity, and the coefficient k is the 
stiffness coefficient. Physical examples of systems with negative stiffness can be 
found in [TT|rT^. 

Equation (jl.l[) is also written in matrix form as 

:')'=^(^)(:')' -lit))- ^'-'^ 

Note that the roots of the auxiliary equation of (|l.ip and the eigenvalues of the 
matrix A are A± = ( — 6 ± \J\P- — 4A;) /2. For constant coefficients 6, k with fc < 0, 
one eigenvalue is negative and one is positive. This makes the point m — 0, u' = 
a saddle point and the zero solution unstable (see definition below). 

The literature for this differential equations with time-varying coefficient has 
several results about instability with negative stiffness, but none about stability. In 
an attempt to extend the stability results to time- varying coefficients, the so called 
frozen coefficient method has been developed. In this technique the coefficients 
are frozen in time and the system is analyzed as a system of constant coefficients 
[U Sec. 10.7]. Thus arises a belief that if the eigenvalues corresponding to time- 
varying coefficients change slowly with respect to time, then the instability obtained 
for constant coefficients remains valid. However, we did not find a precise statement 
of how small should be the rate of change of the eigenvalues. In this article, we 
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disprove the believed instability by showing that for each positive number, there 
exist coefficients b and k for which (jl.ip is stable and the rate of change in the 
eigenvalues does not exceed the given number. See Remark 13.61 

The main objective of this article is to find conditions on b{t) and k{t) for sta- 
bility in the negative stiffness case. More precisely, we find conditions for the 
transition from instability (at the frozen state) to stability of systems with variable 
coefficients. To this end we use Lyapunov functionals in sections 2, and a fixed 
point argument in section 3. Also we show that for every stiffness coefficient, there 
is a damping coefficient that makes p.ip stable. Similarly, for every non-negative 
damping coefficient, wc find a stiffness coefficient so that (11.11) is stable if and only 
if an integral condition (I3.1|l on b is satisfied. Then as an application, we extend 
the stability results to nonlinear systems. We conclude this article by presenting 
some instability results that complement those in the literature. 

In this article, we assume that b{t) and k{t) are continuous functions so that 
standard arguments in differential equations guarantee the existence and uniqueness 
of a solution u{t) = u{t, uq, ui). 

Definition. The zero solution is stable if for each e > 0, there exists a corresponding 
5{e, to) > such that u(to)^ + "'(^o)^ < S'^ implies u{t)'^ + u'{tY < <? for all t > to. 
Equivalently, max{|w(io)|, |M'(to)|} < S implies max{|u(i)|, < e. 

The zero solution is asymptotically stable if for some S > 0, the condition u(to)^ + 
-u'(to)^ < S"^ implies limt^oo = 0. Equivalently, max{|M(to)|, |M'(to)|} < 

S implies limt^oo = limj^oo \u'{t)\ = 0. 

The zero solution is strictly stable if it is stable and asymptotically stable. A 
solution that is not stable is called unstable. For linear systems, the stability of 
one solution implies the stability of all solutions, in which case the system is called 
stable. 



Stability for (jl.ip . with positive stiffness, has been established when b{t) and 
k{t) are bounded above and below by positive constants in [10^. Assuming that |&|, 
I A; I, and \k'\ are bounded above, Ignatiev ^ proved uniform asymptotic stability, 
under the assumption that k and k' /{2k) -\- b are bounded below by two positive 
constants. 

Stability for systems of the form x' = y4(t)x and x' = A{t, x)x has been studied 
by several authors O [H [71 [H [9] . However, their restrictions on the matrix A{t) do 
not allow for negative stiffness. In [T, the eigenvalues have negative real part, and 
the matrix satisfies some growth conditions. In 9J the average of the real part of 
the eigenvalues is negative and the matrix satisfies some growth conditions. 

Our first stability result reads as follows. 

Theorem 2.1. The zero solution of (|1.1|) is stable if, for all t > to, the following 
conditions are satisfied: 

b{t) > 0, —r — + k{t) > M for a constant M (which maybe negative)^ (2.1) 



2. Stability using Lyapunov functionals 




fc)V(26). 



(2.2) 



Proof. First, we define the Lyapunov functional 



E{t) 
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and compute its derivative along the solutions of (jl.ip . 

E'{t) = 4-[e^+'=]u2 + 2(e^+'= - k)uu' - 2h{u'f . 
at 

Then, factoring —26 in the last two terms, and completing the square, we have 

E'it) = [^(e^+'= - kf + j^ei+'y - HYb^ei^' ~ ^> (2.3) 

By (|2.2p . the coefficient of v? is non-positive, and because h> 0, E'{t) < 0; there- 
fore, E{t) < E{to) for aU t > to- 

Now, we show that the zero solution is stable. For e > 0, we select 6 > such 
that 

(5^ < min{e*^ l}e^/ max{l, exp(— ^ + k{to))} ■ 
Note that u{to)'^ + u'(taY < S"^ implies 

min{e^^ l}e^ > max{l, exp(--i- + k{to))}6'^ > E{to) . 
Also note that 

Eito) > E{t) > e^^^uitf + u'{tf > min{e*^ l}{u{t)^ + u'{t)^) . 

The stability of the zero solution follows from the two inequalities above. This 
completes the proof. □ 

Uniform stability is obtained under the additional assumption that -t- fc is 
bounded above; because, the delta in the proof can be chosen independent of to- 

Remark 2.2. The above theorem makes the transition from instability (at the 
frozen state) to stability take place. However, when k{t) is non-positive and non- 
decreasing. Conditions (|2.ip and (|2.2p imply b{t) growing exponentially, which is 
very restrictive. To prove this remark, note that for non-decreasing fc, we have 
(!/&)' < {1/b + ky. Then by (gT]), e*^ < eH'^' and -k> e^^ Then by ((O) . 
(l/fe)'e*^ < -556^^. This implies b' > 6/2, which in turn implies b{t) > 6(0)e*/2. 

Example. Among the equations with exponential damping, there are stable equa- 
tions that satisfy and some others that do not satisfy (|2.ip - (l2.2p . For instance, if 
k — ~1 and 6 = e"*, then (jl.ip has solutions of the form 

exp (^^) (c, [!(-- + -, z) + -, z)] + [/C(- --, z) - /C(- + z)] 

where ci and C2 are the integration constants, z = e"*/27i, and I(a, z) and B{a, z) 
are the modified Bessel functions of first and second kind, respectively. When n = 1 
this solution converges uniformly to ci, when n > 1 it converges to 1c\^pripK ^ as 
i ^ 00. In both cases this implies stability. The zero solution is also stable for 
n < 1/2. For instance, for n = 1/2 the solution converges to 32ci, and for n = 1/4 
it does to 781250ci. However, (|2.ip - (l2.2p are satisfied only for n > 2. 

In an attempt to weaken the growth restrictions on 6, we consider now the case 
where the damping is a positive constant. Note that the assumptions below restrict 
the stiffness to remain negative. Also note that the larger (smaller) the constant 
damping is, the slower (faster) the negative stiffness is needed to be compensated. 
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Theorem 2.3. The zero solution of (II. 1|) is stable if bit) is a positive constant and 
for all t > to: 

There exists a positive constant a such that — k(t) > a, (2.4) 
-k' + 2k^/b<0,. (2.5) 
Proof. Wc define the Lyapunov functional 

E{t) = ~k{t)u{tf +u'it)^, 
whose derivative along solutions of (jl.ip is 

E'{t) = -k'u^ - 2kuu' + 2u'u" 
= -k'u^ - ikuu' - 2b{u'y 
fc2 



= l-k' + 2^]m2 _ 2b[-u + uf. 
b b 

By (|2.5p . the coefheient of is non-negative. Since 5 > 0, E'{t) < so that 
E{t) < E{to) for all t > tg. To show stability of the zero solution, for each positive 
e, we select 6 > such that 

max{l, —k{to)}S'^ < min{l, a}e^ . 

With this delta, we can show that the definition of stability is satisfied, and hence 
the proof is complete. □ 

Note that (|2.5p implies k' > 2k'^/b which yields a lower bound for the rate of 
change in k. Since b is constant, this inequality provides bounds for the rate of 
change in the eigenvalues of A{t): 

— i < , < and — — > — , > . 

dt - b^b"^ - Ak ~ dt - by/b"^ - 4fc ~ 

So that the transition to stability happens when both eigenvalues approach zero 
sufficiently fast. 

3. Stability using a fixed point theorem 

In this section we eliminate the restriction that b must grow exponentially, by 
using a fixed point argument similar to those used in [2]. We start by stating a 
condition that is necessary (but not sufficient) for stability. 

Lemma 3.1. Assume k{t) < for t > tg. Then the condition 

e -''o ''ds < oo (3.1) 

to 

is necessary for stability of (jl.ip . 

Proof. First, using the integrating factor exp ( J^*^ b{s) ds) we transform p.ip into 
the equivalent equation 

(e^'^u')' + ke^''u^O (3.2) 
which is used for setting up a contrapositive argument. Let the initial values u{to) 
and u'{to) be positive. Then by the continuity of the solution there is a non-empty 
maximal interval [ioj^i] where u(t) > 0. On this interval, 

(e/V)' = -ke^^ > 0. 
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So that e-^ ^u' is non-decreasing; hence, 

u'{t) > u'(to)e"^'o'' > 0. 

Therefore, u{t) is increasing and the maximal interval can be extended to [to,oo). 
Integration on the above inequality yields 

u{t) > u{to) + u'ito) / •'^'0 ''ds, yt>to. 

J to 

Note that when (|3.ip is not satisfied, the solution u{t) is unbounded which implies 
being unstable. This completes the proof. □ 

Remark 3.2. For each damping coefficient b{t), there is a stiffness coefficient k so 
that (jl.ip is stable if and only if (|3.ip is satisfied. In fact, 

= - exp ( - 2 / b{s)ds) 

J to 

leads to exp(± J^^e ^^o'' ds) being solutions of (|I.ip . To check stability, we use 

that all solutions can be written as u{t) — cie^ + C2e~'^ with r{t) ~ J^*^ e -^'o ds, 
and that by (|3.ip . r(t) and r'(t) are bounded. On the other hand if (|3.ip is not 
satisfied, one of the two solutions is unbounded which leads to instability. 

To check that the two functions above are solutions, let (j){t) = e'^*-\ Then 
= r'e" and = [r'fe" +r"e\ Since r' = e-/^ {r'f = e-^/^ = -fc, and 
r" = -be-f\ it follows that kcj) + [r'fe'' = and bcf)' + r"e'' = 0. Therefore, is 
a solution of (jl.ip . 

As an illustration of the remark above, we have the following two examples: 
Firstly, when b{t) = 2/(^ + 1), (|3.ip is satisfied and the equation 

u"{t) + j^^u'it) - ^ 

has solutions of the form u{t) — cie^/'-*^^-' + C2e^^/(*+^^ Since both exponential 
functions are bounded on [0,cxd), we can prove stability. Secondly, when b[t) = 
l/{t-\- 1), (|3.ip is not satisfied and the equation 

has solutions of the form u{t) = ci{t + 1) + C2/{t + 1). Since the first function is 
unbounded on [0,cxd), we have instability. 

Now we set up a mapping whose fixed points are solutions of (jl.ip . From (|3.2p . 
it follows that the solution u{t) satisfies 

u'{t) = u'(<o)e" - ^'0 [ k{T)e^'l ''u(t) dr , 

J to 

and 



/ e ■' 


'*o ^ lis — / e 




'to 


J to 


J to 



u{t)^u(to)+u'{tQ) I e J*o"ds- / e J*o " / k{T)e^'o "u{t) dr ds := F[u]it) 

(3.3) 
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For e > 0, we define the following convex subset of the space of continuous differ- 
entiable functions. Let 

/■oo 

B, = {u: \u{tQ)\ < e/4, \u'{to)\ < min{e/4,e/(4 / e"^*o''dr)}, 



\u{t)\ <e, \u'{t)\ <eyt>ta}. 

Using the suprcnium norm ||m||oo = supoj^^ can show that this set is 

closed under the norm HuHoc + ||w'||oo- Note that under assumption p.ip this set 
is not empty; at least, there are constant functions in this set. 
For the next result we define the hypotheses: 

'■'''o''<2 Vs>to, (3.4) 



e 



■'''o'' / \k{T)\e^'^o''dTds < (3.5) 



to "'to ^ 

e^^*o / |fc(T)|e-''*o ''dT < - Vs > ■ (3.6) 
Jto 2 

Lemma 3.3. Under assumptions p.ip and (|3.4p - (j3.6p . f/ie transformation F maps 
into -Be and has a fixed point. 



Proof. Let it be a function in B^. Then by (|3.4p and (|3.6p . 

|F[M]'(t)| < |w'(io)|e"^'o'' + ||7.||^e"^4'' / |A:(T)|e^'o^dT<e/2 + e/2 = e. 
By dsn, 

|FM(t)| < \u{t^)\ + \u'{t^)\ f\-^'o' + \\u\\^ f'e-^^l" r \k{T)\J^o' dT ds 

< e/4 + e/4 + e/2 = e. 

Therefore, F maps into itself. To find a fixed point for F, we define an iter- 
ative process that can start at any function wq in B^. For n — 1,2,..., define 
Wn = F[wn^i]. Note that the values w„(io) and w'„{to) remain unchanged in these 
iterations. Also note that by (|3.5p and (|3.6p . 

e^-^to '' / |A:(T)|e'^'o '' \\wn - Wn-i\\ oo dr ds < -||w„-w„_i||oo 



and 



'Jn)i't)\ < e'-^*o / |A;(T)|e'^'o ''\\wn - w„_i|loo dr < -||w„ - w„-i|| 



Therefore, is a contraction and converges to a fixed point of F in B^; hence 
the solution of p.2p is in the set B^^ and satisfies the conditions for stability. This 
completes the proof. □ 

We are ready to present the main result of this section. 

Theorem 3.4. Under assumptions (|3.ip and p.4p - p.6p . the zero solution of (|l.ip 
is stable. 

Proof. For each e > 0, we define B^ as above and set 5 = min{e/4, e/ (4 e~ ^ '')}. 
For initial conditions |u(to)| < 5 and |u'(to)| < 5, the solution is obtained as a fixed 
point of F; therefore, \u{t)\ < e and \u'{t)\ < e which impUes stabiUty of (|l.ip . □ 
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Remark 3.5. For each stiffness coefficient k{t), there exists a damping coefficient 
b{t) that makes stable. In fact, for a constant a > 1, we let 

b{t) > 2\k(t)\{t + + 2/{t + a) , (3.7) 

so that the conditions in Theorem 13.41 are satisfied with to = 0. Condition (|3.4|) 
follows from b{t) > 0. Note that 6 > 2/(t + a) and J* b> J^,* 2/(< + a) = 2 ln(< + a). 
So that 

exp(- b) < exp(-2 \n{t + la)) = (i + a)^^. 
Jo 

Condition p.ip follows from integrating in the inequality above. Note that from 
((3J)) . <{t + a)~^b/2 - (t + a)'^ which implies 

\k\e^'' < ^{t + a)-^be^ - {t + a)-'^e^ 
Integrating on [0,t], we have 

J^\k\ef'dr<l{t + a)-'ef'-^. 

Then 

JO 2 2 

which is (|3.6p . Integrating on [0,oo), we have 

which is (|3.5p . Therefore, is stable with this choice of b{t). 

Remark 3.6. The believed instability is disproved as follows: For each e > 0, we 
find b{t) and k{t) such that (II. ip is stable and the rate of change in the eigenvalues 
of A(t) is less than e, in absolute value. 

Let b{t) = 4:/{t + a) and k{t) = + a)^ , where a = max{l, ^/bje}. Note 

that p.7p is satisfied, and hence the conditions for Theorem 13.41 are satisfied; so 
that (|l.ip is stable. 

The rate of change in the eigenvalues of A{t) is 

dX± 1 / b \db 1 dk 

~dr^2\ y'b^-ik^di^ y'b^-4^k'dt' 

The absolute value of the coefficient of db/dt is bounded by 1, while \db/dt\ — 
4/{t + a)"^ < i/a^. The coefficient of dk/dt is bounded as follows 

I 1 I _ I 1 ,1/2 ^ 1 _ ft + g) 

Since dk/dt = 3/{t + a)*, and a > 1, 

,dA±, 4 3 4 3 19 
— -\ < 1 < \ < < e 

' dt ' - a2 4(f + a)3 - a2 ^ 4a3 " 

Which proves the claim of this remark. 

We conclude this section with a stability result for the non-linear case. 
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Remark 3.7. As an applications to non-linear equations, we consider the differ- 
ential equation 

u"{t) = f{t,u',u), (3.8) 

where /(t, 0,0) = and / is differentiable at {t,0,0). Then u{t) = is a solution 
of (|3.8p . The stability of the zero solution is studied by considering the linearized 
version 

u"{t) = f2{t,0,0)u' + Mt,0,0)u, (3.9) 
where f2(t,x,y) — dxf{t,x,y) and f^(t,x,y) = dyf{t,x,y). Note that this linear 
approximation is valid only for small values of u and of u' . The stability of (|3.9p is 
studied by setting 

6(t) = -/2(t,0,0), fc(i) = -/3(t,0,0) 

and applying results from this section, without any further modifications. How- 
ever, the instability results in the next section may not hold because the linear 
approximation is valid only for values u, u' close to zero. 

4. Instability using Lyapunov and Chetaev functionals 

Instability of (|l.ip . with negative stiffness, was obtained by Ignatiev [5], assuming 
that \k'\ are bounded above, and that —k and |^-|r + &| are bounded below by 

two positive constants. In the same article, instability is proved when jb^ + fc < 0, 
and when ^iP' + k > Q with some additional assumptions. 

Our next result states that the zero solution is unstable, for non-positive stiffness 
and non-positive damping. 

Theorem 4.1. If b{t) < and k{t) < for all t > to, then the zero solution is 
unstable. 

Proof. For each pair of initial values u{to) > and u'{to) > 0, by continuity of the 
solution and its derivative, there exists an interval where u{t) > and u'{t) > 0. 
Since 6 < and fc < 0, from p.ip . it follows that u"{t) > 0; i.e., u is concave 
up and u'{t) > u'{tQ) > on this interval. Let ti be the largest value such that 
u"{t) > for to < t < ti. If <i < +00, the graph of u is concave up on [io,^i]i 
u{ti) > u{to) > and u'(ti) > u'{to) > 0. By continuity of the solution and its 
derivative, there exists t2 > ii, such that u{t) > and u'{t) > on [ti,t2]. Since 
6 < and k < 0, from (jl.ip . it follows that u"{t) > on [ti,i2]- This contradicts 
ti being maximal; therefore, ti — +oo. 

Because ti — +oo, the graph of u is concave up and u'{t) > u'{to) > for all 
t > to- Therefore, limt^ao u{t) = +oo for all arbitrarily small and positive initial 
values. This implies instability of the zero solution and completes the proof. □ 

Example. The equation 

u"{t) - ^u'{t) - ^u{t) = (4.1) 

satisfies the conditions of Theorem 14.11 and has solutions of the form 

u(t) ^ citX(2, 2Vi) + C2tJC{2, 2Vi) . 
For large t the dominant terms in each branch of this solution are 

Cl A.2Vt C2^ ^ 2^t 
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which diverges as t cx3, which imphes instability. 

Theorem 4.2. The zero solution of (jl.ip is unstable under the following condi- 
tions: For all t > to, 

there exists a positive constant a such that — k{t) > a, (4-2) 
k'{t) + 2b(t)k{t)>0. (4.3) 

Proof. We construct a functional similar to the one in Chetaev's theorem \6[ . How- 
ever, the proof for the variable coefhcient case is not the same as the constant case. 
Let 

vit)^u{tr+^^u'itf, 

whose derivative along the solutions of p.ip is 

V'{t) = 2uu' - ^{ur + T^'u" = -^^±^{u'f . (4.4) 
fc^ k k^ 

Note that V{t) can be positive or negative and that we can select u(to) and M'(io) 
so that V{to) < 0. Then by (|43l) . V'{t) < and V{t) is non-increasing, which 
allows only two possible cases: 

Case 1. When \imt^QoV{t) = —oo, since V > -^(w')^ — ~q("')^i have 
Ihnt^ociu')'^ — +00. Therefore, u' is unbounded and the zero solution is unstable. 
Case 2. When V{t) is bounded below, being non-increasing, it converges to some 
negative number —L^. Then there exists ti such that V{t) < — L^/4 for t > ti. By 

(1121), 

--{uT < tJttK)' < V{t) < -iV4 for t > t,, 
a kit) 

which implies |'u'(t)| > Ly/a/2 > 0. If u' is positive, then it is bounded below by 
a positive constant for t > ti. Therefore, limt_,oo u{t) — +oo and the zero solution 
is unstable. If u' is negative, then it is bounded above by a negative constant for 
t > ti. Therefore, Imit— ,00 u(t) = —00 and the zero solution is unstable. This 
completes the proof. □ 

Theorem 4.3. The zero solution of (jl.ip is unstable under the following condi- 
tions: For all t > to, 

there exists a positive constant a such that — k{t) > a, (4-5) 
k'{t) + 2b{t)k{t) <0, (4.6) 
there exist positive constants and tj such that 

1 (4.7) 

-k'{t) - 2b{t)k{t) > -a^b{tfk{tf, fort>t3. 

Proof. We define the Chetaev's functional and compute its derivative as in the 
proof Theorem 14.21 Next we select u{to) and u'{to) so that V{to) > 0. Then by 
(j4.6p . V'{t) > and V{t) is non-decreasing, which allows only two possible cases: 

Case 1. When limt_^oo Vit) = +00, since V{t) < u{t)'^, we have limi_,oo u{t) — 
-\-oo and the zero solution is unstable. 

Case 2. When V{t) is bounded above, being nondecreasing, it converges to some 
positive number i^. Therefore, Vit) dt converges. Since the integral 1/tdt 
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diverges, a limit comparison yields limt_,oo tV'{t) — 0. From ()4.4p . (|4.5p . and (14. 7|) . 
we have 

= lim tV'{t) > lim a36(t)2(u'(i))2 > 0. 

t — ^oo t — *oo 

Then limt^oo \bu'\ = 0; therefore, there exists a time ti such that \b{t)u'{t)\ < La/4: 
for all t > ti. 

Since V{t) < u{t)'^ and limt^ooV{t) — L^, there exists a time t2 such that 
\u{t)\ > 3L/A for all t > t2- There are two possible cases for t > max{ti,t2}: 
Case 1: u{t) > 3L/4. From and ([431l . 

„ , , , , , La k'dL La 3aL La 

u"{t) = -hu' -ku> ; — > + — — = — > 0. 

4 4 4 4 2 

The solution u, being bounded below and having concavity greater than a posi- 
tive constant, must have limt^oo 'u(i) = +oo. This implies instability of the zero 
solution. 

Case 2: u{t) < -3L/4. From (HHI) and g^), 

,, , . , , , La k3L La 3aL La 

u"{t) = -hu' -ku< \ < = < 0. 

^ ' -4 4-4 4 2 

The solution u, being bounded above and having concavity less than a negative 
constant, must have limt^oo u{t) = —00. This implies instability of the zero solu- 
tion. 

Since the above reasoning applies to arbitrarily small positive initial conditions 
such that ~k{to)uito)^ > u'{to)^, the zero solution is unstable and the proof is 
complete. □ 
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